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TRANSIENT  SKIN  EFFECTS  IK  EXPLODING  WIRE  CIRCUITS 


ABSTRACT 

The  transient  response  of  the  coaxial^  current-measuring  shunt 
commonly  used  in  high-current,  high  frequency  applications  (up  to  1  me) 
is  analyzed  by  Laplace  transform  methods.  An  approximate  solution  is 
obtained  which  allows  estimates  to  be  made  of  the  errors  expected.  The 
current  shunt  measuring  a  damped  oscillation  will  always  report  an  initial 
current  slope  of  zero,  and  the  maximum  rate  of  current  rise  is  sensed 
shortly  after  sv/itch-on.  It  is  several  percent  low  in  typical  cases. 

At  the  first  current  maximum,  the  shunt  reading  is  a  few  tenths  percent 
high  and  lags  z':ie  impressed  current  by  a  small  fraction  of  a  cycle.  The 
transient  resistance  of  an  idealized  plate  condenser  is  analyzed  using 
the  asymptotic  solution  for  current.  A  numerical  calculation  indicates 
no  alteration  of  initial  conditions  on  the  damped  oscillation  to  arise 
from  this  source,  so  long  as  the  characteristic  damping  time  of  the 
transient  skin  effect  is  small  compared  with  the  ringing  time. 
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1.  RESPONSE  OF  A  COAXIAL  CURRENT  MEASURING  RESISTANCE 


The  current  pulse  encountered  in  typical  exploding  vire  experiments 

has  a  transient^  irregular  wave  form  with  its  basic  frequency  in  or  near 

1|. 

the  megacycle  range  and  peak  amplitude  frequently  greater  than  10  amps. 
One  desires  to  measure  both  current  magnitude  and  pulse  shape  with  high 
accuracy  (within  2$  of  peak  values)  by  means  of  a  device  which  itself 
should  add  negligible  impedance  to  the  circuit  and  negligible  pertur¬ 
bation  to  the  current  pulse. 

Such  a  device  woiild  seem  to  be  found  in  the  coaxial  shunt  design 

given  by  Park^^^^  which  has  found  use  in  measurement  of  large  currents 

at  power  frequencies  and  more  recently  in  high  current  generators  at 

(2 ) 

about  60  kc''  '  and  in  exploding  wire  experiments  up  to  20  kc  basic 
frequency. 

Current  shunts  of  this  design  are  presently  in  use  for  the 
measurement  of  current  pulses  thru  exploding  wires  in  a  circuit  with 
ringing  frequency  of  kk6  kc.^^^  While  it  might  appear  that  the  shunt 
technique  is  sufficiently  well  understood  that  no  difficulties  would 
occur  in  extending  the  method  to  l/2  mc^  certain  discrepancies  have 
been  observed  in  the  measurement  of  damped  oscillations  in  the  exploding 
wire  circuit.  These  deviations  from  expected  values  have  pointed  to 
the  necessity  for  better  understanding  of  shunt  characteristics  under 
transient  conditions,  and  have  thus  led  directly  to  the  present 

m 

theoretical  study. 

In  the  modified  design  considered  here  (see  Figure  l),  the  current 
measuring  resistor  consists  of  a  short,  hollow  cylinder  of  nichrome  V 
made  the  inner  conductor  at  the  shorted  end  of  a  coaxial  stub  which 
itself  is  connected  into  the  circuit  in  series.  Voltage  sensing  is 
accomplished  by  a  small  co-ax  line  entering  through  the  shorted  end, 
inside  and  coaxial  with  the  nichrome  element.  In  a  particular  example 

-5 

employed,  the  resistance  of  the  nichrome  element  is  about  10  ohm  and 
its  calculated  inductive  reactance  less  than  %  of  the  resistance  at 


l/2  EC.  The  nichrome  tube  is  thin  enough  to  indicate  nearly  constant 
resistance  independent  of  skin  effect  at  frequencies  up  to  1.5-  me  yet 
strong  enough  to  resist  magnetic  pinch  forces  from  currents  up  to  10^ 


When  a  study  is  made  of  the  damped  oscillations  vhich  can  occur 

vith  the  exploding  wire  replaced  by  a  shorting  bar,  certain  characteristic  ■ 

discrepancies  are  noted.  The  damping  constant  a  =  R/2L  can  be  determined 

from  comparison  of  successive  .current  peaks  and  the  angular  frequency 
2  2  2  2 

(0=  0)  -  a  where  oi  =  1/LC,  from  the  time  intervals  between  successive 

0  0 

zeroes  of  the  current  curve.  If  capacity  C  is  determined  by  an 
independent  bridge  measurement,  the  RLC  constants  of  the  circuits  can 
be  obtained  without  any  recourse  to  actual  current  magnitudes.  Yet 
when  the  values  so  obtained  are  used  to  predict  the  damped  oscillation, 
one  finds  that  the  observed  values  of  initial  current  slope  are 
invariably  smaller  than  the  values  predicted  from  circuit  theory  by 
amounts  several  times  larger  tha.n  the  expected  experimental  error. 

The  theoretical  treatment  of  the  shunt,  skin-effect  problem,  by 

Silebee  and  Park^^^  is  based  on  the  assimption  of  steady-state,  continuous 

wave  operation.  Because  of  this  limitation  no  conclusion  can  be  drawn 

from  their  calculations  as  to  the  accuracy  of  a  coaxial  resistor  for 

meas’uring  switch-on  transient  pulses  like  those  encountered  in  the 

(?) 

e>:ploaing  wire  phenomenon.  On  the  contrary,  the  study  by  Haines  '  of 
inverse  skin  effect  for  a  solid  cylindrical  conductor  carrying  a  steady 
oscillation,  shows  Lmportant  differences  between  the  initial  current 
d:.stribution  and  that  vnhich  exists  after  transients  have  damped  out. 

Thus  the  possibility  exists  that  the  response  of  a  coaxial  current- 
mea^’uring  shunt  may  be  seriously  in  error  at  switch-on  because  of 
tr..nsient  skin  effect.  This  point  of  view  is  given  some  additional 


iort  by  the  work  of  haninger 


who  studies  the  current  distribution 


sccurring  in  fine  wires  when  applied  current  is  a  step  function  rising 
L'rom  -  A  to  A  at  a  given  time. 


0 


l/2  Kc.  The  nichrome  tube  is  thin  enough  to  indicate  nearly  constant 

resistance  independent  of  skin  effect  at  frequencies  ,ip  to  1.5  me  yet 
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strong  enough  to  resist  magnetic  pinch  forces  from  currents  up  to  10 
amps . 

When  a  study  is  made  of  the  damped  oscillations  vhich  can  occur 

with  the  exploding  wire  replaced  by  a  shorting  bar,  certain  characteristic ■ 

discrepancies  are  noted.  The  damping  constant  a  =  R/2L  can  be  determined 

from  comparison  of  successive  current  peaks  and  the  angular  frequency 

oi  =  m  -a.  where  o)  =  1/LC,  from  the  time  intervals  between  successive 
0  o' 

zeroes  of  the  current  curve.  If  capacity  C  is  determined  by  an 
independent  bridge  measurement,  the  RLC  constants  of  the  circuits  can 
be  obtained  without  any  recourse  to  actual  current  magnitudes.  Yet 
when  the  values  so  obtained  are  used  to  predict  the  damped  oscillation, 
one  finds  that  the  observed  values  of  initial  current  slope  are 
invariably  smaller  than  the  values  predicted  from  circuit  theory  by 
amounts  several  times  larger  than  the  expected  experimental  error. 


The  nheoretical  treatment  of  the  shunt,  skin-effect  problem,  by 

Silsbee  and  Park'"^'  is  based  on  tne  assumption  of  steady-state,  continuous 

wave  operation.  Because  of  this  limitation  no  conclusion  can  be  drawn 

from  their  calculations  as  to  the  accuracy  of  a  coascial  resistor  for 

meas’oring  switch-on  transient  pulses  like  those  encountered  in  the 

(?) 

exploaing  v/ire  phenomenon.  On  the  contrary,  the  study  by  Haines  ' ^  of 
inverse  skin  effect  for  a  solid  cylindrical  conductor  carrying  a  steady 
oscillation,  shoves  iirportant  differences  between  the  initial  current 
distribution  and  that  which  exists  ai'ter  transients  have  damped  out. 

Thus  the  possibility  exists  that  the  response  of  a  coaxial  current¬ 
measuring  shunt  may  be  seriously  in  error  at  switch-on  because  of 
tr..nsient  skin  effect.  This  point  of  view  is  given  some  additional 
^u'.port  by  the  work  of  Maninger  who  studies  the  current  distribution 
occurring  in  fine  vires  when  applied  current  is  a  step  function  rising 
from  -  A  to  A  at  a  given  time. 


In  the  present  paper  we  consider  the  transient  response  of  a 
coaxial  resistive  element  to  a  current  suddenly  applied.  Our  treatment 
is  closely  modeled  after  that  of  Haines  who  employs  the  Laplace 
transform  to  obtain  a  general  solution  for  the  case  of  a  solid 
cylindrical  wire  under  arbitrary  current  excitation.  While  it  appears 
that  the  tubular  shunt  problem  can  be  exactly  solved  in  principle, 

certain  difficulties  arise  in  the  analysis  which  make  it  expedient  to 

\ 

employ  an  approximate  solution  to  obtain  usable  results  in  the  practical 
cases  considered.  When  this  is  done  estimates  can  then  be  given  of  the 
accuracy  to  be  expected  from  the  coaxial,  current  measuring  resistance 
during  the  early  portion  of  a  transient  current  pulse. 
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MALYSIS 


2.1  Formulation  of  the  Problem 


Consider  the  diffusion  of  magnetic  field  into  a  long  cylindrical 
tube  of  metal  of  inner  radius  r  and  outer  radius  a.  We  eii5)loy 
cylindrical  polar  coordinates  with  z  taken  along,  and  r  perpendicular 
to  the  axis  of  the  tube.  Then  following  Haines'  treatment  we  assume 
Maxwell's  equations  (mks  units),  neglect  displacement  currents,  set 
B  =  pH  and  use  Ohm's  law,  i^  =  oE^,  in  a  form  which  ignores  both 
inertial  effects  of  the  current  carriers  and  thermal  diffusion.  The 
relevant  current  equation  for  the  i^  is  then 


where  t  is  time,  p.  =  permeability  and  a  the  conductivity 

of  the  medium.  If  v;e  assume  i^  to  be  independent  of  0  and  z,  which 
is  reasonable  on  account  of  symmetry  and  because  we  wish  currents 
anywhere  along  the  tube  to  be  the  same,  then  conservation  of  charge 
plus  the  boundary  conditions  imply  radial  current  everywhere  zero. 

For  convenience  introduce  the  change  of  variables 
2 

X  =  r/a  and  y  =  t/paa  .  Then  Fq.  (l)  becomes 


(2) 


which  is  a  form  of  the  diffusion  equation.  The  transformation  on  r 
is  linear;  whereas  the  time  axis  is  compressed  or  expanded  non-linearly 
depending  on  the  values  of  p,  a  and  a.  Since  the  fundamental  solutions 
will  be  expressed  in  terms  of  x  and  y,  it  is  clear  that  to  reach  a  given 
condition  of  current  penetration  will  take  much  longer  in  media  of  large 
conductivity  and  dimension  than  when  these  values  are  small.  Put 
another  way,  the  solutions  we  shall  obtain  are  self-similar  in  terms 
of  the  scaling  parameter  paa  .  We  require  solutions  to  Eq.  (2)  under, 
the  follov/ing  boundary  and  initial  conditions. 
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with  solution 


i(x,p)  =  A(p) 


I  (p^/^x) 
0 


+  B(p)  K^(p^/^x) 


(8) 


where  I  and  K  are  modified  Bessel  functions  of  the  first  and  second 

®  °  (9) 

kind  foi  imaginary  argument.  The  functions  A(p)  and  B(p)  can  be 
obtained  from  the  boundary  conditions  (5)  and  (6)  after  one  takes  the 
Laplace  transform  of  these  equations  under  the  assumption  that  the 
order  of  operations  can  be  interchanged.  The  result  is 


fn  1 

I(x,p)  =(X 

\d  (q)) 

2 

2na 

1/2  > 


1/2^ 


l/2^ 


K^(ep'^^'^)l  I^(xp'^^^)  xdx  -  I^(ep'^''‘^)  I  K^(xp^/^)  xdx 


(9) 


where  primes  indicate  differentiation  with  respect  to  the  argument. 
The  desired  current  function,  i^(x,y),  can  now  be  found  using  the 
inversion  integral  theorem.  Keeping  Eq,  (9)  in  view  we  write 


i.(^;y) 


a 


1 

iiri 


Lim 

CO 


a 


-ip 


i(x,p)  dp 


(10) 


where  a  lies  to  the  right  of  all  of  the  singularities  of  i(x,p).  The  task 
of  exhibiting  the  function  i^(x,y),  once  the  Laplace  transform  of  the 
applied  total  current  (y)  has  been  determined,  is  now  seen  to  be  that 
of  evaluating  the  open  contour  integral  (lO).  If  the  singularities  of 
the  integrand  are  simple  poles,  then  a  well  known  expansion  theorem  allows 
e;q)ression  of  i^(x,y)  as  the  sum  of  the  residues  of  the  integrand.  This 
is  the  form  in  which  Haines  gives  the  exact  solution  to  the  skin  effect 
problem  for  a  solid  cylindrical  conductor.  No  equivalent  exact  solution 
has  been  found  for  the  tubular  conductor  represented  by  Eq.  (lO).  A 
partial  investigation  of  the  poles  of  Eq.  (9)  shows  that  there  are  none 
in  the  iimiiedlate  vicinity  of  the  origin.  With  the  aid  of  the  definitions 


© 


« 
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and  recurrence  relations  among  the  Bessel  functions  the  denominator  of 

Eq^.  (9)  can  be  transformed,  into  a  more  familiar  form  (see  Eq.  (28)  below) 

whose  first  few  zeros  are  tabulated  in  Jahnke  and  Emde.^^*^^  These 

identify  the  first  members  of  a  denumerable  infinity  of  poles  of  the 

integrand  on  the  negative  real  axis.  Additional  values  c^  be 

calculated  from  the  appropriate  asymptotic  expansions  for  the  Bessel 

functions.  While  the  values  of  the  real,  negative  poles  can  be 

determined  to  sufficient  accuracy,  convenient  analytical  expressions 

for  the  residues  have  not  been  obtained;  furthermore,  no  pro.of  is  at 

hand  excluding  the  possibility  of  complex  poles  of  the  denominator  in 

the  left  half  plane.  Because  of  these  lacks  the  exact  solution  for 

i  (x,y)  cannot  be  written  down  in  explicit  form  convenient  for 
z 

calculation.  While  it  is  reasonable  to  suppose  that  these  gaps  in 
the  analysis  can  eventually  be  closed,  the  difficulties  ere  formidable 
and  make  a  more  accessible,  approximate  method  seem  preferable  for  the 
present.  Accordingly  we  proceed  now  with  the  development  of  such  a 
method, 

2.2  Approximate  Solution 

An  appropriate  simplification  can  be  found  by  examining  Eq,.  (7)  and 
recalling  that  in  order  to  minimize  skin  effect  the  tubular  shunt  must 
be  thin.  This  means  that  the  variation  of  x(=  r/a)  is  typically  limited 
to  the  range  O.9  ire  ^1.0,  and  in  the  case  with  which  we  are 
principally  concerned  the  ratio  e  of  inner  to  outer  radius  of  the  shunt 
is  0.95-  Ideally  the  measuring  element  of  the  shimt  should  be  as  thin 
as  possible,  but  the  necessity  to  resist  the  self  pinch  forces  of  the 
current  requires  greater  wall  thickness  than  is  strictly  needed  for 
electrical  reasons  alone. 

With  this  limitation  on  x  we  see  that  the  coefficient  l/x  in  Eq.  (7) 

/  * 

could  reasonably  be  replaced  by  1/x^  where  g  /^x^  ^  1.  Once  this  is 

*  The  possibility  of  this  simplification  was  pointed  out  to  the  author 
by  Mr.  J.  E.  Hunsicker  of  this  laboratory. 
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done  the  approximate  solution  for  transformed  current  is 


i(x,p)  =  A^(p)  e^l^  +  B^(p)  , 


where 


fi\"  ■  'sr 

Uj 


1  + .  /  1 


and  X  is  taken  to  be  (l  +  £)/2,  Using  the  boundary  conditions  to 
in 

evaluate  and  we  obtain 


i(x,p)  =  cLlu'(q). 

.  2 


-  e)  _  f^e^2^^  “ 


f^e^l(^  -e)  _ 


While  Eq.  (ij)  is  analytically  much  more  transparent  than  its 
counterpart  Eq.  (9);  evaluating  the  inversion  integral  is  still  far 
from  simple.  We  carry  the  approximation  procedure  two  steps  further. 

First;  we  replace  %  in  the  integrals  of  Eq.  (15)  by  the  mean  value 
x^  =  (1  +  €)/2;  and  second,  after  integration,  approximate  the  resulting 
exponentials  by  their  series  expansions,  retaining  terms  up  to  fifth 
degree  in  the  arguments  to  ensure  sufficient  accuracy.  The  choice  for 
x^  may  be  justified  by  noting  from  Eq.  (3)  that  when  current  is  assumed 
constant  over  the  cross  section,  the  area  S  is  given  by  product  of 
mean  circumference  with  thickness  in  the  form  S  =  2rta  [j[l  +  e)/2j 


a(l  -  e)  , 

-  • 


from  which  we  identify  x  as  the  dimensionless  mean 

m 


radius. 


After  the  integration  we  have 


2rta  X 


/l(x  -  e) 


f^e^2^ 


1)  -  ' 


X-.)] 

-  if 


and  upon  expanding  the  exponentials,  keeping  Equation  (l2)  in  mind, 
i(x,p) 


1  1 

r  ; 

[p(x  -  e)^/2| 

[1  -  a  -  £)/3x^....^ 

2rta^x  (1  - 

m 

e)  ^1+1 

Lp(i 

-  .m 

1  -  (1  - 

e)Ax^+  (1  -  e)2/20x/..]j 

(15) 

Notice  here  that  the  denominator  of  the  first  factor  is  just  the  area  S 
of  the  tube  cross  section;  and  that,  aside  from  those  introduced  by  the 
transform  of  the  impressed  current  C^(y),  there  is  only  one  pole;  and  it 
is  located  on  the  negative  real  axis.  It  is  completely  determined  by  the 
thickness  and  mean  radius  of  the  tube. 

2.3  Current  Response  of  the  Tubular  Shunt 

We  can  further  specialize  Eq.  (15)  by  setting  x  =  e,  thereby 
obtaining  only  the  transformed  current  density  at  the  inner  radius 
where  the  sensing  probe  is  located.  This  is  all  we  need  for  comparison 
of  the  shunt  response  with  the  arbitrary  impressed  current. 

Suppose  that  impressed  current  is  taken  to  be  a  damped  sine  wave  of 

unit  amplitude,  l.e.,  cl(y)  =  e”  sin  my.  Then  =  m/[(p+Ci:)^-HO^ 

and  the  Laplace  transform  of  introduces  two  poles  into  Eq.  (I5) 

in  addition  to  the  one  already  found.  When  the  inverse  transform 

is  taken,  three  terms  will  be  obtained,  one  corresponding  to  each  pole, 

and  two  will  be  steady  state  terms  from  the  impressed  current.  The 

third  will  be  a  transient  term  arising  from  the  denominator  of  Eq.  (l5). 

Whereas  the  poles,  and  corresponding  transient  terms  of  the  exact 

solution  Eq.  (9),  are  at  least  denumerably  infinite;  the  approximations 

by  which  Eq.  (15)  is  obtained  reduce  this  infinity  to  one.  Accordingly 

only  one  transient  term  is  to  be  expected  in  the  expression  for  true 

current;  however,  this  transient  may  be  regarded  as  fairly  typical  of 

the  infinity  of  terms  in  the  exact  solution.  Note  that  because  of  our 

I  2  2 

choice  of  the  variable  y  it  follov/s  that  a  =  a  pba  and  m  =  m'  (laa 

I  I  22  1/2 

where  a  =  B./2L  and  o)  =  (i/LC  -  R  /hL  )  '  in  terms  of  the  known 

circuit  constants. 
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Under  the  above  conditions  Eq.  (15)  is  now  simplified  enough 
so  that  an  inverse  transformation  can  easily  be  found.  It  is 

i  (e  y)  =  _ _ 

S  [(1  -  ahf  . 


e”  +  (7/00)  (1  -  Oi/y)  e“  ^  sin 


my 


-  oy 

-  e  ''  cos  my 


d  ' 


(16) 


where 


-1 


7  = 


6/(1  -  e)  J 


1  -  (1  -  e)/h\  +  (1  -  e)  /20x 


m 


At  the  origin  of  time  y  =  0  and  i  (e,0)  =  0. 

The  quantity  S  i^(ejy)  is  the  total  apparent  current  as  reported 
by  the  shunt  on  the  assumption  that  the  current  density  at  x  =  e  is 
constant  over  the  entire  cross  section.  This  apparent  current  is  just 
the  quantity  we  wish  to  compare  with  the  impressed  current  <^(y). 

Certain  general  conclusions  can  be  reached  by  further  analysis  of 
Eq.  (16). 

If  the  shunt  thickness  approaches  zero,  then  e— ^1  and  the  pole  -  7 

takes  on  a  negative  value  as  large  as  we  please.  The  first  and  third 

terms  of  (16)  are  negligible  and  v/e  recover  exactly  the  impressed 
-  (TV 

current  e  sin  my.  Therefore  the  approximate  treatment  predicts 
that  a  sufficiently  thin  shunt  will  measure  the  impressed  current  as 
accurately  as  desired. 

Under  slightly  less  restrictive  conditions,  if  m<<.7  and  1  -  a/y 
is  of  order  unity,  the  first  and  third  terms  have  coefficients  much 
smaller  than  the  dominating  central  term.  Near  current  peak  sin  my~l 
and  the  first  and  third  terms  may  be  neglected.  The  amplitude 
of  the  sin  my  term  is,  to  good  approximation,  just  1  +  (7/7;  so  an 
estimate  of  accuracy  of  the  peak  current  reading  can  be  made.  In 
the  examples  of  interest  it  is  also  true  that  a <<7  and  magnitude 
of  expected  error  in  the  current  peak  is  less  than  1^. 


Ih 


2.k  Inll^al  Current  Slope 

To  examine  the  initial  current  rise  sensed  by  the  shunt  we  need 
expressions  for  the  first  and  second  derivatives  of  Eq.  (l6).  With 
the  definition  I  =  S  i  (e,y),  these  may  he  written  as 


dl 

£ 


tse 


-  oy 


(1  -  a/rf  +  (Wyf  L 


and 


d^I 


me 


-  ay 


2  2 
(i  -  a/r)  +  (Wy) 


cos  my  -  e“  +  (m//  -  a/m  +  c^/oyy)  sin 

(17) 


(m/7  -  2a/m  +  a  /7m)  m  cos  my  +  70  ''' 

2  2  ^2 

-  (1  +  a/7  -  a  /m  ■>r  cc lt£>  j)  sin  my  . 

(18) 


From  Eq..  (17)  it  is  clear  that  when  y  =  0  the  first  derivative 
vanishes;  so,  regardless  of  thickness,  the  shunt  reports  initial  slope 
as  zero,  in  contrast  to  the  positive  value  m  for  the  impressed  current. 
While  the  slope  of  the  impressed  cxirrent  falls  monotonically  to  zero  in 
the  first  quarter- cycle,  the  slope  measured  by  the  shunt  must  rise  to 
some  maximum  value  and  then  decrease.  This  maximum  is  necessarily  less 
than  m,  although  for  the  class  of  thin  shunts,  in  which  we  are  mainly 
interested,  the  difference  is  only  a  few  percent.  The  argument  proceeds 
as  follows. 


The  maximum  of  slope  will  occur  where  the  second  derivative,  Eq.  (I8), 
vanishes;  thus  the  bracket  on  the  right  must  be  zero.  Since  this  must 
occur  while  my  l/k,  i.e.,  quite  early  in  the  cycle  if  the  shunt  is  to 
reproduce  the  wave  form  at  all  accurately,  we  can  approximate  the 
trigonometric  terms  by  1  and  my  respectively.  If  a<m  <'<7,  the 
condition  reduces  to 

(m/7  -  2a/m)  m/7  -  (m/7)^  7y  +  e“  ^  =  0  ,  (19) 
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from  vhich  an  approximate  solution  can  be  readily  found.  Given  the 
solution  as  y  =  y^,  Eq.  (19)  and  the  inequalities  can  be  used  to 
simplify  the  expression  for  the  maximum  slope  sensed  by  the  shunt. 
The  result  is 


(20) 


and  we  see  that^  to  this  approximation^  the  indicated  maximum  slope 
is  always  less  than  that  of  the  impressed  current. 

2. ^  Typical  Cases 

In  Table  I  values  are  given  for  two  nichrome,  shunts  used  for 
current  measurements  in  the  recent  series  of  exploding  wire  experiments. 

The  sensing  element  of  the  shunt  in  each  case  had  an  outer  radius 

“3  6 

a  =  3-175  X  10  mj  conductivity  a  =  O.926  x  10°  mho/m^ 


TABLE  I 

6  a  0)  7  7y^  1  +  a/7  l-(a/7)(7y3_+2) 

Shunt  I  .926  3-ii95  32.8it  1117.  k.66k  I.003  0.979 

Shunt  II  .880  Same  Same  L30.1  1-008  O.956 

2 

and  the  time  scaling  factor  |j.cra  =  11. 72  psec.  The  last  two  columns 
show  the  deviations  in  peak  current  and  initial  slope,  respectively, 
to  be  expected  on  account  of  the  transient  skin  effect. 

Measurements  of  initial  current  slopes  from  experimental  records 
provided  by  the  shunts  may  yield  values  less  by  several  percent  than 
the  estimates  given  in  the  last  column  of  Table  I.  In  the  practical 
case,  an  average  is  taken,  graphically  or  analytically,  among  measured 
slopes  near  the  origin  of  -Dime.  Because  of  the  characteristic 
response  of  the  tubular  shunt  given  by  Eq. (16)  r  (18),  values  of  slope 
less  than  the  maximum  will  predominate  and  the  experimentally  determined 
value  must  necessarily  be  less  than  the  true  maximum  the  shunt  can 
indicate.  Slope  discrepancies  up  to  105^  may  conceivably  be  accounted 
for  in  this  way. 
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If  the  experimenter  weights  the  limiting  slope  value  at  the 
apparent  origin  in  time  most  heavily^  under  the  impression  that  the 
maximum  slope  will  appear  when  t  =  y  =  0,  then  a  value  may  he  obtained 
even  more  in  error  than  if  he  gives  equal  weight  to  all  early  slopes. 
Failure  of  the  oscilloscope  to  follow  accurately  the  rise  of  shunt 
voltage  will  introduce  an  additional  factor  tending  to  produce  too 
low  an  initial  slope. 

In  making  the  preceding  estimates  it  is  presumed  that  the  scaling 
2 

factor  iioa  is  accurately  known,  but  an  unsuspected  error  in  its  value 
would  have  a  proportionate  effect  in  the  values  of  a  and  cu  without 
changing  the  value  of  7.  While  radius  a  is  easily  determined  to 
within  1^,  finding  the  appropriate  values  of  ^  and  a  for  high, 
cirrrent-rate  pulses  presents  much  greater  difficulty.  Increases 
in  a  and  o)  corresponding  to  an  increase  in  n,  say,  will  tend  to  reduce 
the  value  of  but  not  as  rapidly  as  the  ratio  a/y  increases,  and 
increase  both  the  deviation  in  peak  current  and  the  lowering  of  initial 
slope  from  the  expected  values. 

In  concluding  this  section  we  repeat  again  that  transient  skin 
effects  will  necessarily  cause  a  tubular  current  shunt  to  report  low 
values  of  Initial  current  rise  and  the  error  may  be  10^  or  larger. 
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TRAI^SIEMT  RESISTMCE  OF  CIRCUIT  ELEMENTS 


From  the  nature  of  the  diffusion  of  an  electromagnetic  wave  into  a 
conductor^  it  is  intuitively  clear  that  the  initial  resistance  of  such 
a  conductor  must  be  infinite  at  the  time  of  incidence  of  the  wave. 
Thereafter  the  resistance  falls  rapidly  to  the  value  characteristic  of 
the  ordinary  skin  effect  in  the  case  of  a  continuous  oscillation  or  to 
the  d.c.  value  if  no  oscillation  is  involved.  The  question  arises  under 
what  conditions j  if  any,  can  this  transient,  excess  resistance  seriously 
alter  the  initial  conditions  determining  a  d.anped  oscillation.  Simple 
calctilations  with  the  damped  sine  wave  formula  indicate  that  if  the  circuit 
resistance  is  larger  than  its  steady  state  value  by  a  factor  of  2.5 
in  the  first  quarter-cycle,  the  first  current  peak  will  be  too  low  by 
25^0.  It  follows  that  higher  average  resistance  over  a  shorter  time 
interval  can  accomplish  the  same  resilLt.  Interest  then  focuses  on 
the  exact  form  of  time  dependence  of  the  transient  resistance;  for  a 
resistance  varying  hyperbolically  with  time  as  current  increases 
linearly,  can  prodnce  a  finite  RI  drop  au  t  =  0  and  thus  alter  the 
initial  conditions  under  which  a  damped  oscillation  would  commence  and 
in  particular  reduce  the  initial  ciu-rent  slope. 

Since  most  of  the  resistance  of  a  typical  discharge  circuit  is 
concentrated  in  the  condenser  bank,  we  shall  cxairane  the  resistance 
variation  of  an  ideal  flat  plate  condenser  vath  respect  to  transient 
skin  effect  and  provide  estimates  ol^the  change  in  initial  conditions 
which  m.ay  occur  from  this  cause.  The  trea.tment  will  consist  in 
solution  of  the  flat  plate  condenser  problem  and  use  of  the  asymptotic 
solution,  valid  near  the  origin  in  time,  to  compute  numerically  the 
resistance  variation  of  a  condenser  similar  in  some  respects  to  those 
actually  employed.  It  will  also  be  shown  that  the  asymptotic 
solutions  for  the  solid-cylindrical  and  cy] ir.rii’ical-tube  conductors 
are  identicrd  in  fuictional  form  with  that  for  the  flat  plate;  so 
that  similar  behavior  caun  be  expected  for  connecting  wires,  spark 
terminals,  etc.  'which  have  cylindrical  symmetry.  On  account  of  the 
co.mparatively  large  dimensions  of  these  elements  their  contribution 
to  the  total  resistance  is  small  and  may  be  considered  negligible. 


3.1  Flat  Plate  Condenser  -  Asymptotic  Solution 


We  seek  the  ciirrent  distribution  in  an  infinite  flat  plate  of 
thickness  a  representing  one-half  of  a  parallel  plate  condenser. 

I 

Assume  rectangular  coordinates  vith  x  in  the  direction  transverse 

1 

to  the  plate  and  z  in  the  direction  of  cxurent  flow.  As  before, 

X  and  y  are  dimensionless  space  and  time  variables  such  that 
X  =  X  /a  and  y  =  t/iina  with  0  £•  x  £:  1.  The  entire  formulation 
given  in  §  2.1  may  be  paraphrased  and  the  analysis  carried  through 

t 

in  detail.  The  analog  of  Eq.  (5)  in  terms  of  i^(x,y)  (where  z  =  z  ) 
can  be  written  asj^i^/Sj^  =  0.  In  order  for  the  total  current, 

to  remain  finite  it  is  necessary  to  consider  a  finite  portion  of  the 
plate  of  width  w.  One  finds  without  difficulty  that 


i(x,p) 


_  l/p  cosh  Cvp(l  -  x)3  X'f q)j' 
wa  sinh  /p 


(21) 


We  investigate  the  asymptotic  representation  of  this  solution,  and 
its  inverse  transform,  for  large  p  corresponding  to  small  values  of 
y.  If  p  is  sufficiently  large  the  Laplace  transforms  of  any  of  the 
following:  current  ramp,  sine  wave  or  damped  sine  wave  are  all  the 
same,  viz.. 


I 


(22) 


and  we  take  the  condition  for  the  validity  of  (22)  to  be  p 
an  inequality  which  has  aiready  arisen  in  connection  with  the  approximate 
solution.  If  (1  -  x)  >  2.5  and  to  avoid  an  infinite  p  we  take  x  <  1, 
i.e,,  X  is  bounded  away  from  1,  the  hyperbolic  functions  in  Eq.  (2l) 
can  be  replaced  to  an  accuracy  of  about  lj>  by  their  leading  exponential 
terms .  Thus 


I(x,p) 


“L 

wa 


(23)  . 
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with 


/p  7  2. 5/(1  -  x) 
P  7?  0)  ?  q:  . 


(23a) 

(23b) 


and 

While  in  the  analysis  up  to  Eq.  (23),  p  is  considered  real  and 

positive,  when  the  inversion  is  considered  via  contour  integral 

Eq.  (10),  the  analytic  continuation  of  Eq.  (23)  throughout  the 

entire  complex  plane  is  involved.  A  branch  point  is  seen  to  exist 

at  p  =  0  and  the  inverse  is  found  by  integrating  along  both  sides  of 

—  ^ 
a  contour  which  when  shrunk  coincides  with  the  negative  real  axis. 

To  estimate  the  upper  limit  on  y^O  for  which  the  asymptotic  solution 
might  be  valid,  we  observe  that  the  exponential  in  Eq.  (lO)  will  damp 
out  any  contribution  to  the  integral  occurring  along  the  negative  real 
axis-  save  in  the  range  0  ^  |ypj  ill-  The  upper  limit  of  order  unity 
is  somewhat  arbitrary,  but  would  probably  never  be  chosen  larger  than 
5.  Thus  we  expect  the  asymptotic  solution  to  be  valid  for  j  p  |  large 
and  Oi:  y<  l/|  p  [ .  The  Inequalities  {2%)  and  (23b),  or  equivalent  ex¬ 
pressions,  serve  to  define  the  smallest  permitted  magnitude  of  p  and 
therefore  the  largest  permitted  value  of  y. ■ 

With  these  conditions  in  mind  we  now  quote  the  inverse  of  Eq.  (23) 
from  a  table  of  transforms.  The  inverse  is 

^  (V™)  [2  exp  (-  x^/4y)  -  X  erfc  (x/2  l/y)  ,  (24) 

where 


erfc  (x)  =  1  -  erf  (x)  =  (2/ 


*  Cf .  the  discussion  given  by  Carslav/  and  Jaeger,  Conduction  of  Heat 
in  Solids  2nd  Ed.,  (Oxford  University  Press,  1959)  PP  302-303  and 
Appendix  I. 
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An  alternate  form  is 


=  Wwa) 


(2l;a) 


We  note  In  paaeing  that  if  x  ^  0,  i^  and  all  of  its  derivatives  vith 
respect  to  y  vanish  at  y  =  0  because  of  the  singularity  of  the  function 
at  this  point. 


3.2  Resistance  Variation  with  Time 

To  find  the  time  dependent  resistance  of  the  flat  plate  we 

calculate  the  electrical  power  dissipated  in  a  section  of  length  I 

and  width  w.  With  the  help  of  Ohm's  law  and  the  current  density 

i  we  find 
z 

1 

P(y)  =  (iwa/ff)  j  (x^y)  dx  .  (25) 


We  may  also  express  the  power  as  product  of  Instantaneous  resistance 
times  total  current  squared.  Thus 


P(y)  =  R(y)  I  (y) 


-T  2 


=  R(y) 


wa 


(x,y)  dx 


(26) 


Equating  Eq.  (25)  and  (26)  and  using  =  i/owa  leads  to 


r  / 

-  1 

R(y)/R^c  "  \  ^ / 

\  (x,y)  dx 

-^0  / 

Jo 

(27) 


an  expression  that  is  free  of  the  dimensions  of  the  plate.  For 
this  reason  it  can  be  taken  as  representative  of  the  resistance 
variation  of  the  entire  condenser. 
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Using  Eqs.  (21+)  and  (27),  numerical  calculations  have  been  made 

for  an  aluminm;  flat  plate  condenser  with  plate  thickness 

a  =  1.1+  X  10  m,  the  same  as  that  of  the  foil  in  the  condensers 

7  ,  2  -9 

actually  employed.  For  aluminum,  a  =  3*7  ^  10  mho/m  and  poa  =  9*1  x  10 

(6) 

sec.  In  the  exploding  wire  studies  previously  mentioned^  '  the  basic 
frequency  is  0.1+46  me  and  the  appropriate  value  of  o)  is  .0255.  If 
we  adopt  this  value  of  tu  as  representative,  it  is  clear  that  any 
value  of  p  larger  than  unity  will  satisfy  inequality  (23b)  quite  well. 

On  the  other  hand  (23a)  imposes  a  much  more  severe  restriction  on 
values  of  p  and  y.  For  example  if  x6’0.5  the  requirement  (23a)  demands 
P/'25  and  the  relevant  interval  of  y  is  0— y^.04  which  corresponds  to  a 
time  interval  of  less  than  l/2  nano-sec.  It  might  be  thought  that 
integrating  over  the  range  O^xjrl  in  Eq.  (27)  would  involve  such  an 
inaccurate  representation  of  i^  in  the  neighborhood  of  x  =  1  that  the 
calculation  would  be  vitiated.  Such  is  not  the  casej  because,  for  these 
early  times  (y^.04),  the  current  is  concentrated  mainly  in  the  region 
xi:0.5  and  the  contribution  to  the  integrals  from  the  remainder  of  the 
range  is  typically  a  few  tenths  of  a  percent  and  in  the  worst  case  less 
than  The  calculations  show  that,  in  this  range  of  y  values 

2.8^R(y)/R^^  i:23.4,  i.e.,  transient  resistance  falls  steeply  from 
infinity  to  about  three  times  the  dc  value  in  l/2  nano-sec.  In  this 
interval  the  resistance  variation  is  represented  within  about  5^  by 
the  relation  R(y)/K^j,  =  1  +  0.377/y^'^^^* 

The  significance  of  this  result  lies  in  the  implication  that  for 
a  current  increase  linear  with  time  and  consequently  linear  with  y, 
the  limit  as  y-^0  of  R(y)  i  (y)  is  zero.  Thus  the  Ri  term  in  the 
circuit  equation  for  a  damped  oscillation  has  no  limiting  value  other 
than  zero  and  cannot  influence  the  initial  rise  of  current.  For  this 
reason,  we  can  eliminate  the  transient  condenser  resistance  as  a 
possible  cause  of  a  slope  discrepancy  at  t  =  0.  Furthermore,  the 
fall  from  infinite  resistance  to  values  comparable  with  the  steady  state 
resistance  occurs  in  a  small  fraction  of  the  first  quarter- cycle.  This 
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calculation  shows  R(y)/R^^<l.l  by  t  =  ^.1  nano-sec.  which  interval 
is  only  2<j,  of  the  first  quarter-cycle.  Therefore^  while  the 
calculation  for  y/'.Oli-  is  admittedly  of  problematic  accuracy  because 
the  limiting  inequalities  (23a)  and  (25b)  are  no  longer  satisfied, 
the  indications  are  strong  that  transient  resistance  cannot  contribute  't 
any  appreciable  deviations  from  the  expected  damped- sine-wave  behavior 
of  the  circuit. 


3.5  Asymptotic  Solutions  -  Cylindrical  Geometry 

We  complete  our  discussion  by  indicating  here  the  form  of  the 
asymptotic  solutions  characteristic  of  the  tubular  current  shunt  and  of 
the  solid  cylindrical  conductor.  We  outline  the  procedure  omitting 
details . 

Starting  from  Eq.  (9)  and  making  use  of  differential  relations 

(5) 

among  the  I's  and  K's,''^'  one  obtains 

"  2na^  [Ii(  Vp)  K“(r/p)  -  I^(e  v/p)  K^( '/p)J  • 

(28) 

Then  with  the  asymptotic  expressions  for  the  Bessel  functions  with 
arguments  greater  than  10  (Cf.  ref.  5  PP-  220-221),  using  Eq.  (22) 
and  dropping  all  but  the  leading  terms,  Eq.  (28)  simplifies,  with 
an  accuracy  better  than  %,  to 

I(X,P) - g - cosh  >^(x  -  ,  (29) 

27fa  l/x  p  '  sinh  V^(l  -  e) 

under  the  restrictions  that 

\/p7  2.5/(1  -  e)  (29a) 

and  p^lCtD^a  .  (29b) 
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If  the  inequality  /^(x  -  e)  >2.5  can  be  maintained,  then  we  can  write 


- 2 - 

2na  l/x 


x) 


(30) 


and  the  close  analogy  with  Eqs.  (23),  (23a)  and  (23b)  is  immediately 
seen.  As  before,  an  infinity  of  p  will  occur  to  preserve  the  inequality 
in  the  neighborhood  of  x  =  e;  but,  as  in  the  previous  case,  the 
transformed  current  distribution  will  be  least  in  the  immediate  vicinity 
of  the  critical  value  x  ='  e,  with  the  probable  result  that  integrals 
for  current  and  power  will  be  only  slightly  affected  by  the  inaccurate 
representation  offered  by  (30). 

In  a  similar  way  we  can  show  without  difficulty  that  the  asymptotic 

(3) 

form  for  the  cylindrical  conductor  obtained  from  Haines'  Eq.  (ll),'  '  is 

exactly  Eq.  (30)  where  inequality  j/p>12.5  (l/x  +3)  must  hold  if  1^ 

accuracy  is  desired.  For  this  case  x>0  is  required  to  preserve  the 

inequality,  but  the  usual  arguments  may  be  given  for  using  the 

transformed  asymptotic  current  in  calculations  of  current  and  power. 

In  cylindrical  coordinates  these  integrals  will  contain  an  additional 

factor  of  X,  but  the  behavior  of  the  integrand  will  be  dominated  by  the 

function  given  in  Eq.  (24).  We  may  then  state  with  confidence  that  the 

transient  skin-effect  resistance  in  all  three  cases,  viz.,  planar, 

cylindrical  and  tubular  will  follow  similar  lavs.  We  may  draw  the 

conclusion  that  no  adverse  effect  on  initial  conditions  in  oscillatory 

exploding  wire  circuits  is  to  be  expected  from  transient  skin  effect 

in  condenser,  shunt  or  lead  wires.  These  conclusions  only  hold  when 

2 

the  characteristic  time,  y.aa  ,  is  negligible  compared  with  the  time  of 
one  oscillation.  Should  these  times  become  the  same  within  an  order  of 
magnitude,  then  we  may  expect  appreciable  and  important  effects  to  arise 
from  the  diosipation  caused  by  the  transient  skin  effect  resistance. 


F.  D.  BEKNETT 
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